第 40 卷 第 10 期 
2023410 H 


基于 ANUSPLIN 的 黄土 丘陵 区 气象 要 素 
插值 方法 对 比 与 评估 
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摘 要 : 气象 要 素 是 反映 地 球 水 热 过 程 的 关键 因子 ,气象 数据 的 准确 获取 对 生态 保护 和 农业 生产 研究 具有 重要 
影响 气象 要 素 的 插值 结果 ,制约 了 数据 的 准确 性 。 为 了 获取 黄 
丘陵 区 的 温度 、 降 水 的 空间 分 布 ,探索 地 形变 化 对 气象 空间 插值 结果 的 影响 ,本 研究 基于 专业 气象 插值 软 


义 。 黄 土 丘陵 区 是 典型 的 丘陵 沟 塞 区 ,地 形 严 
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ANUSPLIN 以 延 河 流域 及 周边 2010—2021 年 105 个 气象 站 点 逐日 温度 降雨 资料 为 基础 ,引入 25 m、90 m、1 km 三 个 


不 同 精度 的 数字 高 程 模型 作为 协 变量 进行 插值 ,从 而 获得 黄土 丘陵 区 温度 降水 忆 
区 较 好 的 适 


变化 规律 ,评估 ANUSPLIN 插值 方法 是 否 能 在 黄土 丘陵 


有 格 数据 以 揭示 该 地 区 降水 的 时 空 
j。 结 果 表明 :(1) 延 河流 域 东 部 延长 地 区 温度 


较 高 .西部 温度 较 低 ;降雨 量 空间 插值 显示 中 部 与 西北 部 较 低 ,东部 较 高 ,温度 降雨 丝 符 合 以 往 气象 站 数据 规律 ， 
ANUSPLIN 模 型 对 黄土 丘陵 沟 蜜 区 温度 降雨 空间 插值 有 较 好 的 适应 性 ,(2) 在 三 种 不 同 的 分 辨 率 DEM 模拟 场景 下 ， 
温度 插值 精度 排序 为 :25 m > 90 m > 1 km; 降 雨 插 值 精度 排序 为 :90 m>25 m>1 km。 本 研究 的 结果 为 丘陵 区 气象 分 


布 和 插值 提供 了 参考 依据 。 


关键 词 : ANUSPLIN 模型 ; 气候 要 素 ; 黄土 丘陵 区 ; 数字 高 程 数 据 ; 评估 


气象 因子 是 影响 农业 生产 、 人 类 活动 以 及 生态 
环境 的 重要 要 素 ,尤其 是 在 干旱 半 干 旱地 区 。 气 象 
要 素 在 农业 和 环境 等 相关 领域 研究 较 多 ,主要 应 用 
于 模型 预测 和 作物 生产 等 方面 "。 关 于 气象 的 研 
究 ,最 初 人 们 依赖 于 气象 站 的 雨量 计 监 测 ,但 由 于 
气象 站 点 数目 较 少 ,空间 不 完全 稚 盖 , 且 时 效 性 较 
差 ,在 完整 捕 换 气候 事件 和 描述 空间 变化 等 方面 能 
力 不 足 。 随 着 空间 处 理 和 分 析 技 术 的 发 展 ,空间 插 
值 的 方法 被 用 来 描述 气象 因子 的 空间 分 布 ,一 定 程 
度 上 实现 了 气象 空间 覆盖 ,提高 了 人 们 对 气象 空间 
变化 的 认识 。 

目前 ,插值 方法 较 多 , 且 在 系统 模拟 和 反 演 领 
域 都 有 不 同 程度 的 应 用 ,常见 的 方法 有 局 部 多 项 式 
法 (Local polynomial interpolation), FF 7& P % YE 
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(Spline)、 克 里 金 插 值 法 (Kariging 和 反 向 距离 权重 法 
(Inverse Distance Weight, 简称 IDW) 等 14。ANUSP- 
LIN 是 由 Hutchinson" 开发 用 于 不 同 空 间 和 时 间 尺 
度 水 文 气象 插值 的 方法 ,其 主要 特点 是 可 包含 线性 
协 变量 例如 地 形变 量 等 ,广泛 应 用 于 高 山 丘 陵 等 地 
形变 化 起 伏 较 大 的 地 区 “1。 此 外 ,ANUSPLIN 具 有 
简单 的 操作 流程 ,不 需要 对 其 参数 进行 单独 的 预先 
校准 ,生成 的 温度 降水 时 间 序 列 曲 面具 有 较 高 的 精 
度 和 可 靠 性 外 。 近 年 来 ,ANUSPLIN 插值 方法 已 被 
国内 外 学 者 广泛 应 用 于 气候 数据 的 插值 ”"。 黄 土 
丘陵 区 位 于 黄土 高 原 中 部 ,是 我 国生 态 恢 复 的 重点 
区 域 之 一 ,该 区 沟 餐 纵横 ,地 形 起 伏 多 变 ,温度 降水 
等 气象 因子 受 地 形 影 响 较 大 ” ,传统 的 插值 方法 多 
通过 气象 站 空间 距离 进行 估算 ,而 黄土 丘陵 区 多 以 
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山地 为 主 ,地 形变 化 足以 影响 空间 温度 和 降水 的 差 
异性 分 布 呈 , 单 从 距离 的 远近 获得 气象 数据 并 不 能 
真实 的 反映 该 地 区 的 气象 变化 ,从 而 无 法 满足 植被 
恢复 和 重建 工作 的 需求 。 因 此 ,针对 地 形 复杂 的 丘 
陵 区 ,探索 包含 不 同 分 辩 率 地 形变 化 的 气象 插值 方 
法 对 于 进一步 完善 气象 空间 数据 和 相关 研究 具有 

本 文采 用 ANUSPLIN 空间 插值 方法 ,对 延 河流 
域 2010 一 2020 年 的 温度 降水 进行 空间 插值 ,同时 选 
取 25 m、90 m、1 km 三 个 不 同 精度 的 数字 高 程 模型 
(Digital Elevation Model , DEM) 作 为 协 变量 ,对 比分 析 
各 个 尺度 温度 降水 插值 结果 的 精度 ,评估 ANUSP- 
LIN 插值 方法 在 黄土 丘陵 沟 餐 区 的 适用 性 ,并 选择 
区 域 温度 降水 最 优 数 据 集 , 补 充 丘 陵 区 气象 插值 精 
度 不 够 的 短 板 ,为 准确 获取 地 形 起 伏 较 大 的 丘陵 区 
气象 空间 插值 数据 提供 可 参考 的 理论 依据 和 方法 
支持 。 


1 数据 与 方法 


11 研究 区 概况 

本 研究 以 典型 黄土 丘陵 区 的 延 河流 域 为 研究 
[x dak (36°23’~37°17'N , 108°45'~110°28'E), Z X7 te 
面积 为 7687 km ,海拔 454~1765 m ,该 流域 从 东南 


向 西北 依次 经 过 延长 、 延 安安 塞 和 志 丹 ,是 典型 的 
丘陵 区 (图 1)。 气 候 全 年 较 干 ,夏季 多 十 ,冬季 降水 
较 少 ,处 于 暧 温带 大 陆 性 季风 半 湿 润 气候 向 温带 半 
干旱 气候 的 过 渡 区 ,全 年 降雨 量 496 mm ,年 平均 气 
温 8.8 CC。 从 东南 向 西北 ,降雨 温度 呈现 出 明显 的 
梯度 变化 特征 。 土 壤 以 黄 绵 土 为 主 , 植被 有 森林 
(刺槐 和 辽东 栎 )、 灌 木 从 和 草地 等 。 

1.2 数据 来 源 与 处 理 

1.2.1 气象 数据 基础 气象 数据 来 自 延 河流 域 及 周 
边 站 点 。 数 据 获取 时 间 为 2010 一 2021 年 ,数据 类 型 
为 日 值 数 据 ,包含 日 平均 温度 .最 高 温度 .最 低温 
度 .日 降水 量 等 气候 要 素 , 本 研究 通过 多 年 (11 a) 的 
年 均 温 度 和 年 均 降水 进行 插值 。 由 于 是 长 时 间 序 列 
数据 ,部 分 数据 存在 异常 或 缺失 等 ,因此 ,数据 经 过 
筛选 .剔除 异常 数据 ,最 终 保留 105 个 气象 站 (图 1)。 
1.2.2 DEM 数 据 延 河 流域 23m 分辨 率 DEFM 是 通 
过 对 日 本 ALOS 卫星 搭载 的 PALSAR 传感器 25 m 
DEM 数据 进行 处 理 得 到 的 。ALOS 25 m DEM 可 从 
NASA 的 ASF DAAC (https://search.asf.alaska.edu/#/) 
免费 下 载 。90 m 和 1 km 分辨 率 DEM 数据 来 自 
CGIAR-CSI SRTM 中 国 区 域 数据 。 

1.3 研究 方法 

1.3.1 ANUSPLIN 模型 基本 原理 ANUSPLIN 是 一 
种 使 用 薄板 平滑 样 条 对 噪声 多 变量 数据 进行 透明 
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图 1 延 河 流域 及 周边 气象 站 点 


Fig. 1 Meteorological stations in the Yanhe River Basin and its surrounding areas 
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oy TAAL TA. WOE PEAR PR ESE BE 
是 利用 一 个 平面 去 拟 合 所 有 的 气象 站 点 ,使 得 该 曲 
面 能 通过 站 点 构成 “ 样 条 ”, 得 到 盘 近 所 有 控制 点 弯 
曲 最 小 的 光滑 曲面 , 即 利 用 M 次 多 项 式 ,对 曲线 进 
行 分 段 拟 合 ,用 局 部 的 薄板 来 绘制 分 段 连续 的 曲 
线 。TPS 是 对 样 条 函数 法 的 曲面 扩展 ,插值 过 程 中 
利用 最 优 的 光滑 参数 实现 通 真 度 和 光滑 度 最 佳 平 
衡 ,保证 了 精度 可 靠 , 同 时 插值 曲面 光滑 连续 "1。 
此 外 ,ANUSPLIN 还 允许 转换 自 变量 和 因 变 量 ,并 允 
许 处 理 具 有 缺失 数据 值 的 数据 集 。 当 对 因 变 量 应 
用 变换 时 ,ANUSPLIN 允许 拟 合 曲面 的 反 向 变换 , 计 
算 相 应 的 标准 误差 ,并 纠正 这 些 变换 引起 的 小 偏 


数据 与 插值 结果 进行 比较 ,以 评估 气象 模拟 的 准确 
性 和 偏差 "。 三 种 分 辩 率 的 插值 精度 比较 适合 使 
用 交叉 验证 法 。 本 研究 插值 的 效果 采用 均 方 根 误 
差 RMSE 来 评估 ,表达 式 如 下 : 


(4) 


式 中 :Z. 和 忆 分 别 代 表 气 候 数据 实测 值 和 插值 结果 ; 
n 为 计算 序列 长 度 。RMSE 的 值 越 小 表示 气象 模拟 
效果 越 好 。 

此 外 ,将 实测 值 与 预测 值 进行 相关 性 分 析 ,用 
确定 性 系数 尼 表 示 一 元 多 项 式 回 归 方程 拟 合 度 的 
高 低 。 尼 取 值 分 布 在 0~1 之 间 , 尼 值 越 大 ,表示 模拟 


差 。 人 研究 发 现 ,这 在 将 表面 拟 合 到 降水 数据 和 其 他 
自然 为 正 或 非 负 的 数据 时 特别 方便 *。 

注 板 样 条 估计 z(x) 由 一 个 适当 光滑 的 函数 计 
算 , 该 函数 最 小 为 : 


> 0 


式 中 :y; 为 观测 i 的 观测 数据 值 ;f(x) 为 被 估计 的 样 条 
函数 ;和 为 一 个 正 数 , 称 为 平滑 参数 ;J,(f) 是 用 m 阶 
偏 导数 定义 的 函数 f 的 粗糙 度 的 度量 Lm 阶 偏 导数 
称 为 SPLINA 输入 中 的 样 条 的 阶 ( 它 也 被 称 为 粗 烽 
BEB) I". 

局 部 薄 盘 光滑 样 条 (Partial thin plate smoothing 
splines ) 的 理论 统计 模型 公式 : 

Z,=f(«,) +b y; +e, (i=1, N) (2) 

式 中 :2 为 位 于 空间 i 点 的 因 变量 ;7 为 送 代 次 数 ;f(t) 
为 估算 关于 x 的 未 知 光滑 函数 ,x 为 独立 变量 ;y; 为 p 
维 独立 协 变 量 ;5 为 y 的 p 维 系数 ;e; 为 随机 误差 ""。 

上 式 中 ,也 数 f 和 系数 5 是 通过 最 小 二 乘 估 计 来 
确定 的 : 


| +AJ,(f) (3) 


i ie 
l 
i=l W; 


式 中 :w 为 权重 的 已 知 局 部 相对 变异 系数 ;万 (及 为 函 
数 flo) 的 粗糙 度 测 度 函 数 ,定义 为 7 的 偏 导 ( 称 为 样 
条 次 数 ,也 叫 粗糙 次 数 ); A 为 正 的 光滑 参数 ,在 数 
据 保 真 度 与 曲面 的 粗糙 度 之 间 起 平衡 作用 。 在 
AUNSPLIN 中 通常 用 广义 交 又 验证 GCV 的 最 小 化 以 
及 最 大 似 然 法 GML 的 最 小 化 来 确定 "1。 

1.3.2 插值 精度 评价 指标 ”本 人 研究 选取 数据 的 15% 
作为 验证 集 ,选取 过 程 采用 随机 抽样 ,将 随机 站 点 


效果 越 好 。 
2 结果 与 分 析 


2.1 插值 结果 

基于 ANUSPLIN 分 别 获取 了 研究 区 2010 一 
2021 年 年 均 温和 累计 降水 的 插值 曲面 见 图 2, 从 图 
中 可 以 看 出 , 相 较 于 1 km 的 插值 ,90 m 与 25 m 分 辨 
率 的 插值 结果 细节 表现 得 更 加 突出 ,能 反映 局 部 地 
形 特征 ,也 有 较 好 的 平滑 度 。 

由 图 2 还 可 以 看 出 ,东部 延长 县 逐日 年 均 温 度 
较 高 ,西部 靖 边 . 志 丹 以 及 中 部 的 安 塞 区 逐日 年 均 
温度 较 低 ,三 种 分 状 率 的 逐日 年 均 温度 整体 呈 “ 东 高 
西 低 "的 趋势 。25 m 分 辩 率 的 温度 插值 范围 为 
7.725~13.329 °C ( E| 2a) , 90 m 分 辨 率 为 7.786~ 
13.354 C( Al 2b), 1 km 分辩 率 为 7.872~13.313 C( Al 
2c)。 

从 降雨 插值 可 以 看 出 , 低 降 雨 的 空间 分 布 主要 
集中 在 延 河流 域 西 部 的 靖 边 地 区 .中 部 的 宝塔 区 以 
及 东部 延长 的 少 部 分 区 域 。25 mm 分辨 率 的 逐日 降 
雨 插值 范围 为 1.062~1.847 mm, 90 mm 分 辩 率 的 搬 值 
范围 为 1.046~1.852 mm, 两 者 最 低 降雨 量 相差 
0.016 mm, 上 限 相 差 0.05 mm ,较为 接近 ,但 1 km 分 
辩 率 的 插值 结果 与 23 m、90 m 的 插值 结果 相差 过 
大 ,下 限 为 1.117 mm, 比 90 m 的 结果 高 0.071 mm, HE 
25 m 的 插值 结果 高 0.055 mm; 上 限 为 1.73 mm , 比 两 
者 都 低 , 与 25 m 的 插值 上 限 相 差 0.117 mm ,整体 插 
值 结果 范围 缩小 。 

从 插值 的 温度 降雨 空间 图 中 可 以 看 出 ,细节 特 
征 比 较 突 出 ,有 是 有 明显 的 “马赛 克 ” 纹 路 ,与 过 去 的 
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图 例 
1000 m 逐 日 年 
LC 


图 例 
1000 m 逐 日 年 均 
降水 量 /mm 


= 
1.114 


图 2 基于 不 同 分 辩 率 DEM 的 延 河流 域 温 度 降雨 分 布 
Fig.2 Temperature and rainfall distribution map of Yanhe River Basin based on DEM of different resolutions 


Cokriging 插值 相 比 ,没有 出 现 较 大 的 斑 块 与 “ 牛 
He” ,说明 ANUSPLIN 有 比较 好 的 平滑 度 , 且 很 好 的 
拟 合 了 地 形 特 征 , 这 与 刘 林 等 的 研究 结论 一 致 。 
此 外 ,插值 结果 显示 延 河流 域 的 中 部 、 西 北部 降雨 
较 低 , 东 部 降雨 量 较 高 ,这 在 刘 春 利 等 ”的 研究 结 
果 中 得 以 验证 。 且 土壤 水 分 变化 也 与 之 一 致 弛 。 
温度 插值 特征 也 符合 以 往 气象 站 公布 的 数据 规 
律 。 延 河流 域 的 西北 部 为 草原 地 带 ,整体 呈 降 水 
低 ,温度 低 ,东部 为 森林 草原 地 带 与 森林 带 ,降水 
高 ,温度 高 ,与 插值 结果 一 致 ,说 明 ANUSPLIN 模型 
能 较 好 的 适应 黄土 高 原 丘陵 沟 蜜 区 。 
2.2 插值 结果 比较 

本 研究 以 25 m、90 m、1 km 三 种 分 辨 率 的 海 拨 
数据 为 协 变量 对 黄土 丘陵 沟 殉 区 的 温度 降雨 进行 搬 


值 , 相 较 于 前 期 研究 中 采用 的 1 kmx1 km 或 500 mx 
500 mm 分 辩 率 ” ,本 研究 在 一 定 程度 上 提高 了 精度 ， 
减 小 了 插值 误差 。 从 图 3 可 以 得 出 ,对 于 不 同 DEM 
分 辨 率 下 温度 插值 结果 ,根据 计算 的 均 方 根 误差 
(RMSE) ,25 m 分 辩 率 的 插值 数据 精度 较 好 ,90 mm 分 
辩 率 插值 数据 精度 其 次 ,1 km 插值 数据 精度 较 差 。 而 
在 确定 性 系数 (RR) 方 面 ,25 m 分 辩 率 的 记 为 0.65, 拟 
合 效果 较 好 ,1 km 分 辨 率 拟 合 效果 其 次 ,90 m 分 辨 率 
拟 合 效 果 稍 差 。 综 上 所 述 , 在 三 种 不 同 精度 的 DEM 
分 辩 率 模拟 情景 下 ,温度 插值 精度 排序 为 :25 m> 
90m> 1 km。 贾 洋 等 六 的 研究 表明 ,DEM 分辨 率 越 
高 ,ANUSPLIN 模型 插值 精度 也 越 高 ,温度 插值 结果 
也 遵循 这 一 规律 ,高 精度 的 地 形 数据 可 以 有 效 反映 
空间 差异 ,也 可 提高 温度 空间 插值 的 精确 度 。 
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图 3 基于 不 同 分 辩 雄 DEM 温 度 降水 观测 值 与 预测 值 分 析 结 果 
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Fig. 3 Analysis results of observed and predicted values of temperature rainfall based on DEM with different resolutions 


降雨 插值 结果 精度 检验 按照 年 均 温 来 检验 , 空 
间 数 据 取 365 d。 对 于 降雨 插值 结果 的 RMSE ,90 m< 
25 m<1 km, E} 90 m 分 状 率 的 插值 精度 最 高 ,1 km 的 
插值 精度 较 低 。 同 时 ,90 m 分 辩 率 的 插值 结果 与 实 
测 结果 拟 合 时 ,90 m 分 辨 率 也 是 三 种 分 辨 率 中 最 大 
的 ,R=0.833, 拟 合 效 果 良 好 ,其 次 是 25 m 分辨 率 ， 
1 km 分 辨 率 的 拟 合 效果 最 差 。 这 与 温度 插值 规律 
并 不 相同 ,笔者 推测 出 现 这 一 结果 的 原因 可 能 是 由 
于 黄土 丘陵 区 山地 较 多 天 ” ,降水 容易 在 沟 地 汇集 ， 
导致 站 点 数据 存在 误差 。 尤 其 是 黄土 丘陵 区 普遍 
坡 长 在 70~200 m 之 间 ,90m 的 空间 精度 与 地 形变 化 
精度 恰好 一 致 ,确定 性 检验 仍 需 后 续 工 作 进一步 验 
证 。 但 是 ,该 结果 也 说 明 在 丘陵 区 设置 气象 站 点 要 
考虑 地 形变 化 。 

高 分 辩 率 数据 集 (包括 遥感 数据 和 长 时 间 序 列 
的 气候 数据 ) 是 目前 黄土 丘陵 区 相关 研究 迫切 需要 
的 。 且 近年 来 干旱 区 气候 变化 较 大 加 ,实时 气象 数 
据 对 于 气候 -植被 研究 和 人 类 活动 具有 重要 价值 。 


为 该 区 域 生态 恢复 和 气候 变化 研究 提供 基础 数据 
zg, 


结论 与 讨论 


本 研究 选取 黄土 高 原 典 型 丘陵 区 延 河 流域 作 
为 研究 区 ,以 研究 区 内 及 周边 105 个 气象 站 2010 一 
2021 年 逐日 气象 数据 为 基础 资料 ,以 1 km、90 m, 
25 m 三 种 分 辨 率 的 数字 高 程 数据 为 协 变量 ,利用 
ANUSPLIN 插值 方法 对 研究 区 不 同 地 形 分 辨 率 气 象 
数据 进行 空间 搬 值 ,通过 气象 站 站 点 交叉 验证 方 
法 ,对 比分 析 各 分 辩 率 数字 高 程 数据 的 插值 精度 水 
平 。 人 研究 发 现 数字 高 程 数 据 的 分 辩 率 影响 气象 空 
间 插 值 的 结果 ,具体 结果 如 下 : 

(1) 温度 插值 显示 , 延 河流 域 东 部 延长 地 区 温 
度 较 高 .西部 温度 较 低 符合 以 往 气 象 站 数据 规律 ”; 
降水 量 空间 插值 显示 中 部 与 西北 部 较 低 , 与 该 区 域 
蒸 散发 时 空 变化 趋势 基本 特征 一 致 汪 ,延安 市 周边 


3 


这 一 气象 数据 和 空间 分 布 研究 对 于 干旱 区 气候 变 
化 研究 非常 宝贵 ,同时 高 精度 分 辩 率 的 数据 集 可 以 


本 研究 降雨 量 较 少 ,其 主要 原因 可 能 是 人 为 干扰 时 
致 气象 站 点 数据 存在 误差 ,以 往 研 究 也 表明 , 延 河 
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可 


40 卷 
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流域 的 延安 市 是 人 类 活动 较为 频繁 的 区 域外 ,在 
赵 美 亮 等 在 西北 地 区 水 文 空间 分 布 研究 中 表明 ， 
土地 利用 改变 会 影响 水 文 要 素 分 布 。 总 的 来 说 ， 
ANUSPLIN 模型 对 黄土 据 陵 沟 密 区 温度 降雨 空间 插 
值 有 较 好 的 适应 性 。 

(2) 在 三 种 不 同 的 DEM 分 辩 率 模拟 场景 下 , 温 
度 插值 精度 排序 为 :25 m >90 m> 1 km ,符合 数字 高 
程 数据 分 辩 率 越 高 ,插值 精度 越 高 六 ;降雨 插值 烛 
度 排 序 为 :90 m>25 m>1 km, 不 完全 符合 以 上 结论 ， 
表明 降雨 数据 空间 插值 存在 最 佳 分 辩 率 的 数字 高 程 
数据 ,最 佳 数字 高 程 数 据 与 该 区 域 地 形 特征 有 关 。 在 
野外 观测 中 发 现 , 黄 土 丘陵 区 山脉 坡 长 多 60~300 m 
之 间 ,与 该 结果 具有 一 致 性 ,坡度 和 坡 长 等 地 形 数 
据 可 能 更 适合 描述 降雨 数据 信息 ,而 不 是 25 m 的 海 
PASH o 

本 研究 基于 ANUSPLIN 插值 方法 ,结合 数字 高 
程 数 据 进 行 多 种 空间 尺度 的 空间 精度 验证 ,该 结果 
与 观测 结果 拟 合 较 好 ,说 明了 结果 具有 可 靠 性 ,所 
得 结论 能 够 为 丘陵 区 插值 的 应 用 提供 参考 。 然 而 ， 
本 研究 仅 考虑 了 海拔 数据 作为 协 变量 ,然而 坡 向 、 
坡 位 等 地 形 因 子 对 插值 结果 影响 也 不 可 忽视 ,将 多 
种 协 变量 考虑 在 ANUSPLIN 插值 方法 中 ,可 以 为 气 
象 插值 精度 优化 方法 提供 更 有 意义 的 参考 依据 。 
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using the ANUSPLIN model 


XIAO Xu', ZHENG Cheng’, DING Chengqin’, FAN Chenzhe*, BAI Yuejiang', 
LIN Longchao', YAN Ting', GAO Yu', SHI Haijing™’ 

(1. Yan’ an Meteorological Bureau, Yan’ an 716000, Shaanxi, China; 2. College of Grassland Agriculture, Northwest 
A & F University, Yangling 712100, Shaanxi, China; 3. Institute of Soil and Water Conservation, Northwest A & F 
University, Yangling 712100, Shaanxi, China; 4. College of Natural Resources and Environment, Northwest A & F 

University, Yangling 712100, Shaanxi, China; 5. Institute of Soil and Water Conservation, Chinese Academy of 
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Abstract: Meteorological elements are the key factors used to assess the earth’ s hydrothermal processes. 
Accurate acquisition of meteorological data is thus of great significance to the development of ecological 
protection systems and agricultural research. Loess hilly areas typically have hills and gullies, and consequently, 
the terrain seriously affects the interpolation of meteorological data and its accuracy. Based on the professional 
meteorological interpolation software ANUSPLIN, this study has used the daily temperature and rainfall data 
from 105 meteorological stations in and around the Yanhe River Basin from 2010 to 2021 as the basis, and three 
digital elevation models (DEM) with resolutions of 25 m, 90 m, and 1 km as covariables to interpolate the grid of 
temperature and precipitation in the loess hilly region. The spatiotemporal variation of precipitation was 
analyzed, and the applicability of the ANUSPLIN interpolation method in the loess hilly region evaluated. The 
results show that temperatures in the eastern extension area of the Yanhe River Basin were higher, when 
compared with those in the western area. The rainfall was lower in the central and northwestern regions when 
compared with the eastern regions. The distributions of temperature and rainfall were consistent with the laws of 
the previous meteorological station data. The ANUSPLIN model had a good adaptability to the spatial 
interpolation of temperature and rainfall in the loess hilly- gully region. In three different DEM resolution 
simulation scenarios, the accuracy of temperature interpolation was ranked as 25 m > 90 m > 1 km, and the 
rainfall interpolation accuracy was ranked as 90 m > 25 m > 1 km. 
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